Excitation of stellar oscillations by gravitational waves: 
hydro dynamic model and numerical results for the Sun 
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ABSTRACT 

Starting from a general relativistic framework a hydrodynamic formalism is derived that yields 
the mean-square amplitudes and rms surface velocities of normal modes of non-relativistic stars 
excited by arbitrary gravitational wave (GW) radiation. In particular, stationary GW fields are 
considered and the resulting formulae are evaluated for two general types of GW radiation: radia- 
tion from a particular astrophysical source (e.g., a binary system) and a stochastic background of 
gravitational waves (SBGW). Expected sources and signal strengths for both types of GW radia- 
tion are reviewed and discussed. Numerical results for the Sun show that low-order quadrupolar 
g modes are excited more strongly than p modes by orders of magnitude. Maximal rms surface 
velocities in the case of excitation by astrophysical sources are found to be v < 10 nuns , 
assuming GW strain amplitudes of h < 10 -20 . It is shown that current models for an SBGW 
produced by cosmic strings, with fi G w ~ 10~ 8 — 10~ 5 in the frequency range of solar g modes, are 
able to produce maximal solar g-mode rms surface velocities of 10~ 5 — 10~ 3 mms -1 . This result 
lies close or within the amplitude range of 10~ 3 — 1 mm s _1 expected from excitation by turbulent 
convection, which is currently considered to be responsible for stellar g-mode excitation. It is 
concluded that studying g-mode observations of stars other than the Sun, in which excitation by 
GWs could be even more effective due to different stellar structures, might provide a new method 
to either detect GWs or to deduce a significant direct upper limit on an SBGW at intermediate 
frequencies between the pulsar bound and the bounds from interferometric detectors on Earth. 
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1. Introduction 

Excitation of solar and solar-like oscillations in 
non-relativistic stars is generally attributed to tur- 
bulent convection. Building on first analytical re- 
sults by Goldreich & Keeley ( 1977|, sub sequent 
theoretical studies (e.g., Balmforth||1992[ |Goldre- 



entropy fluctuations. However, the latter term 
was shown to be negligible for excitation of so- 



lar g modes (Belkacem et al. 2009 Appourchaux 



et al.pOlO ). Whereas theoretically computed ex 



ich et al.|[l994 



ISamadi fc Goupiij |200l[ |Belka- 



cem et al. 20081 extended and refined the model 
to describe stochastic excitation of stellar normal 
modes by turbulent convection and identified two 
major driving terms in the wave equation. The 
first source term is related to the Reynolds stress 
tensor and the second term is due to turbulent 



citation rates for radial solar p modes were found 
to be in good agreement with observations by the 
GOLF instrument aboard the Solar and Helio- 
spheric Observatory (SOHO) spacecraft ( |Samadi| 
et al.||2003b|a[ [Belkacem et al.||2006a|b[ ), quantita- 



tive estimates for rms surface velocities of solar g 
modes differ from each other by several orders of 



magnitude (e.g., 


Gough||1985 


Belkacem et al.||2009 


). This is 



to the adopted eddy-time correlation function (see 



1 



Belkacem et al.|2009 and the discussion in |Appour^| 
chaux et al ||2010[ ). 

In this paper, we build on these earlier works on 
excitation mechanisms and on recent work by the 
authors of this paper ( |Siegel fc Rot"hl[2010l hence- 
forth Paper I) and develop a hydrodynamic model 
that describes excitation of normal modes of non- 
relativistic stars by arbitrary, externally generated 
gravitational waves (GWs). In the past, many au- 
thors investigated the possibility of exciting elas- 
tic waves in solid astrophysical objects like planets 



or the Moon by GWs (e.g., Forward et al. 



1961 



Weber|i"967l |Dyson|[i"969l |Zimmerman fc Hellings| 



1980 1 and even excitation of normal modes of fluid 
stars by GWs was considered (|Boughn fc Kuhn] 



1984 Khosroshahi & Sobouti 1997). However, 



many of these models are based on highly sim- 
plifying assumptions. At the advent of asteroseis- 
mology, mainly stimulated by high-precision ob- 
servations of stellar oscillations from satellite mis- 
sions like CoRoT (iBaglin et al.l 12006} and Ke- 



pler ( Christensen-Dalsgaard et al. 2009), we start 



to elaborate on the topic again and present here 
for the first time a thorough hydrodynamic model 
that describes this excitation mechanism properly. 
Our model holds for any non-relativistic star and 
arbitrary external GW fields, it is compatible with 
current models for other excitation mechanisms 
(see above) and it allows us to use realistic cur- 
rent solar and stellar models. In fact, it turns 
out that the formalism presented here is a gen- 
eralization of the formalism by Samadi & Goupil 



(2001) and Belkacem et al. (2008) to include ex- 



citation of normal modes by GWs. Numerical re- 
sults for the Sun that we present in the second part 
of the paper suggest that this excitation mech- 
anism together with experimental data from the 
Kepler satellite and theoretical computations in 
the case of excitation by turbulent convection is 
likely to either provide a new method of detecting 
GWs or of deducing a significant upper limit on a 
stochastic background of gravitational wave radi- 
ation (SBGW) at intermediate frequencies in the 
near future. 

The paper is organized as follows. We start 
from Einstein's field equations of general relativity 
in Section [2] in order to derive an equation of mo- 
tion for the oscillations of a non-relativistic star 
under external forcing by GWs. After present- 
ing the general solution to this excitation prob- 



lem given by an expression for the mean-square 
amplitudes of normal modes, we focus on station- 
ary external GW fields and derive comparatively 
simple analytical formulae for mean-square ampli- 
tudes and intrinsic rms surface velocities in this 
case. In Section [3j these formulae are further eval- 
uated for two general types of possible incoming 
GW radiation: excitation by a particular astro- 
physical source, such as a binary system, and ex- 
citation by an SBGW of astrophysical or cosmo- 
logical origin. In the second part of this paper, 
Section |3J we focus on the Sun as our nearest 
(non-relativistic) star and present numerical re- 
sults for the excitation mechanism developed in 
the preceding sections. In this context, the so- 
lar model and damping rates are discussed. We 
then proceed by reviewing possible astrophysical 
sources of GW radiation, various models for dif- 
ferent candidates of an SBGW as well as expected 
signal strengths. This finally allows us to numeri- 
cally compute mean-square amplitudes and intrin- 
sic rms surface velocities for both types of GW 
radiation; the results of these calculations are pre- 
sented in Sections |4.5.2| and |4.6.2| Although re- 
sults for solar p modes are shown, the focus pri- 
marily lies on solar g modes, and in particular the 
low-order ones, since these modes are shown to be 
more strongly excited than p modes and since p 
modes are by far more readily excited by turbu- 
lent convection. Conclusions are drawn in Section 
[5] and two appendices are added to streamline the 
presentation in Sections [2] and [3j 

2. Hydrodynamic model 

2.1. General relativistic framework 

The starting point for the formation of a hydro- 
dynamic theory that describes excitation of nor- 
mal modes of non-relativistic stars by externally 
generated GWs is, as in Paper I, the general rel- 
ativistic conservation equation of energy and mo- 
mentun^J 

TV. V = 0, (1) 



or equivalently 



(2) 



1 Throughout this paper, Greek indices take spacetime values 
0, 1, 2, 3, whereas Latin indices take spatial values 1, 2, 3 
only. Repeated indices are summed over. 



2 



which follows directly from Einstein's field equa- 
tions of general relativity and the Bianchi iden- 
tities; T^ v denote the components of the energy- 
momentum tensor of the star under consideration. 
Here, we assume an ideal fluid, in which case the 
energy-momentum tensor is given by 



u^u" + pg p 



(3) 



where the symbols p, p, c, and u^ L are used to 
denote density, pressure, the speed of light, and 
the internal four-velocity field, respectively. 

For the reason stated in Paper I, we work in 
Fermi normal coordinates with the origin at the 
center of mass, P(t), at all times. In this reference 
frame, the metric up to second order in the co- 
ordinates x^ is given by the following expressions 



( |Manasse fc Misner||1963||Misner et al.||1973[ ): 



5oo 

9ij 



"I - RolOmX X 
2 RolimX X 



Vij ^RiljmX X 



I rn 



pva p 



(4) 
(5) 

(6) 

de- 



where r)^ v = diag(— 1, +1, +1, +1) and R 
note the components of the Riemann curvature 
tensoiF] evaluated at P(t), that is along the world 
line x 1 = 0. 

In the case of non-relativistic stars, the lin- 
earized theory of general relativity is appropriate 
by definition. Consequently, the metric can be as- 
sumed to be of the form 



9pU ^]pV hfll>i 



(J) 



with h,„, being very small 



! >v Mif., \ v i- v L^iimii, i^^^l 1, and stars 
are taken to move in flat spacetime, i.e., deviations 
from flat space are solely due to GWs, — hr^ , 
and effects of curvature due to the mass of the 
star are neglected. To linear order in h llVl the 
components of the Riemann tensor in linearized 
theory are given by 



Rpvpo — 2 ph per ~\~ ^p^a^up 



- dpdpKa - d v d a h w ). (8) 



2 The conventions for the Riemann tensor, the connection, 
etc., used throughout this paper correspond to the con- 



ventions employed by |Misner et al.| ( |1973} and |Maggiore| 
|2008j. 



It is important to stress that the Riemann tensor 
in linearized theory is invariant, rather than just 
covariant, and thus it can be evaluated in any co- 
ordinate frame. For convenience, we choose h^ v to 
be in transverse-traceless (TT) gauge defined by 



hop = d J h l:j =h\ = 0. 



(9) 



As in Paper I, there are some well-motivated 
approximations to be made. First, we adopt the 
long- wavelength approximation, i.e., it is assumed 
that the typical length scale X/2ir, over which the 
incident gravitational radiation with wavelength A 
changes substantially, is much larger than the di- 
ameter d of the star, d/(A/27r) <C 1. This results 
in h^u having no spatial dependence over the vol- 
ume of the star. Second, we assume that the inter- 
nal motions of the star are purely Newtonian, i.e., 
non-relativistic. Consequently, only terms includ- 
ing T 00 need be retained on the right-hand side of 
Equation with T 00 given by T 00 = pc 2 . 

2.2. Equation of motion 

In order to derive an equation of motion from 
Equation ([!]) or ^ that describes the dynamics 
of the internal velocity field of the star under ex- 
ternal forcing by GWs, the components T M " of 
the energy-momentum tensor and the Christoffel 
symbols have to be evaluated in Fermi normal co- 
ordinates. This is done by expressing these quan- 
tities in terms of the metric components given by 
Equations (|4])-(|6]). Due to their invariance prop- 
erty, the components of the Riemann tensor can 
then be computed with the help of Equations ^ 
and ([9]) . Finally, adopting the aforementioned ap- 
proximations and only retaining terms linear in 
the coordinates x^ and linear in the metric pertur- 
bations hu V yields the usual equation of continuity 
as the zeroth component of Equation ([2| , 



dp 
dt 



v-(H = o, 



(10) 



where v denotes the internal velocity field of the 
star, and a modified Euler equation, 



dpv 
~dt 



V-(pv<S)v) = -Vp + f 



GWi 



(11) 



which originates from the spatial components of 
Equation ^2y, The latter equation has an addi- 
tional term compared to the usual Euler equation 



3 



of hydrodynamics, f GW (x,t), with components 
given by 

/gw = V'/'V"' = Pfhw (12) 

This term represents an additional external driv- 
ing force per unit volume exerted by GWs. De- 



tailed derivations of Equations (10) and (11) are 
given in Appendix |A") 

Effects of the gravitational field of the star it- 
self can now be included by adding a Newtonian 
potential term, pg. Thus, we finally arrive at 

dpv 



dt 



V-(pv (gi v) = — Vp + pg + f 



gw- 



Here, g must satisfy Poisson's equation, 



V-g 



-AnGp : 



(13) 



(14) 



where G denotes the gravitational constant. 



Equation ( 13 ) describes all internal motions of 



a non-relativistic star under external forcing by 
GWs. In the following, we closely follow the pro- 



cedure presented by Samadi & Goupil (2001) in 



order to derive an equation of motion solely for 
the oscillations of the star that are excited in this 
way. As a by-product of applying this procedure, 
we obtain a generalized and unified theoretical de- 
scription of the mode excitation processes due to 
turbulent convection and GWs, combined in a sin- 
gle consistent forced wave equation. However, in 
contrast to 



Samadi & Goupil 



work by Goldreich & Keeley 



(2001) and previous 



(1977) and Goldre- 



ich et al. (1994), which the former is based on, 



we do not adopt the Cowling approximation, that 
is we do not neglect perturbations in the New- 
tonian gravitational potential. Instead, we show 
that terms related to these perturbations can be 
retained as far as oscillations are concerned. 

All physical quantities are decomposed into a 
sum of the equilibrium quantity, q a , and an Eule- 
rian fluctuation, qi, q = qo + q\. We retain terms 
up to quadratic order in the fluctuations pi, p\, s\, 
g 1 , and v. However, since in calculations of stellar 
normal modes only terms up to linear order in the 
perturbations are consistently retained, we neglect 
nonlinear terms that neither appear as damping 
nor as driving terms in the resulting equation of 
motion. Proceeding in this way, the perturbed 
equations that we consider are written as 

dpi 
dt 



V-(pv) = 0, 



(15) 



dpv 
"dt 



Vpi - pig + po9i 

= -V-(p v®v) + pf GW , (16) 

V- 9l = -4irG Pl . (17) 

Additionally, from the assumption of adiabatic os- 
cillations one has 



ds t 
~dt 



dsi 
~dt 



v • V(s + si) 



(18) 



where s t denotes the Lagrangian entropy fluc- 



tuation due to turbulent convection (Samadi & 
GoupilpOl l. 



Equations (15), (17), and (181 are now em- 



ployed to derive the desired equation of motion for 
the oscillation modes from the temporal derivative 



of Equation ( 16 ). The velocity field v is separated 
into a component v osc due to oscillations and a 
component u due to all other motions, the most 
prominent of which is turbulent convection, 



u. 



(19) 



Equation ( 19 1 implies that the turbulent medium 



is decoupled from the oscillation modes. This ap- 
proximation is essential in treating the excitation 
problem analytically, since it is ultimately respon- 
sible for the resulting equation of motion being a 
linear, inhomogeneous and thus analytically solv- 
able differential equation. Furthermore, we as- 
sume incompressible turbulence, V • u = 0. The 
result is the inhomogeneous wave equation for the 
oscillation modes, 



o 2 



l>n \ ^ - £ ) fosc + V(v osc ) 



— Qf^fRey -^cntr + /gw)' (20) 

where the wave operator £ is a linear differential 
operator acting on arbitrary differentiable vector 
fields X according to 

C(X) = — (v [c 2 s V-(p X) + a s X ■ Vs ] 
Po I 



g V-( Po X)\-GV 



V-(poX) 



(21) 

Here, V denotes the volume of the star under con- 
sideration. The linear damping operator T> is iden- 
tical to the one obtained by |Samadi fc Goupfl] 
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(2001) and it is given by Equation (13) therein. 
On the right-hand side of Equation (p0| there are 



three driving terms, the first two of which are, 
respectively, the Reynold source term, which is 
due to turbulent Reynolds stress, and the entropy 
source term, which is due to turbulent entropy 
fluctuations; definitions of these terms are pro- 



vided by Equations (11) and (12) of Samadi & 



Goupil (2001). For the present analysis, however, 



we are solely interested in the third driving term, 



/gw — Pol 



GW' 



(22) 



which represents excitation of normal modes by 
external GWs; here we neglected the contribution 
/°i/gws since it is comparatively small. Conse- 
quently, we simply ignore the first two driving 
terms. This, however, is only possible because of 
the source terms being independent of the terms 
on the left-hand side of the equation, thanks to 
the separation of the velocity field into an oscilla- 
tory and a turbulent component. Hence, the final 
equation of motion to be solved is written as 



-£>(?w) = 



at ' 



(23) 



We note the striking formal similarity of Equation 



( 23 1 with the equation of motion derived in Paper 



I in the case of adopting the isotropic elastic sphere 
model of a star, which corresponds to considering 
the Sun as a spherical resonant mass GW detector. 

2.3. Square amplitudes and surface veloc- 
ities 

This subsection is dedicated to solving the 



equation of motion ( 23 ) and to thereby gaining 



analytical formulae for the square amplitudes and 
surface velocities of global stellar normal modes 
excited by a not yet to be specified external GW 
field, which is represented by the external driving 
force /gw Henceforth, we will drop the subscript 
"osc" and solely use v to denote the velocity field 
due to oscillations. 

We start by discussing the corresponding ho- 
mogeneous equation without driving and damping 
terms, 

d 2 \ 

^-£)v hom = 0, (24) 



in which case the velocity field can be assumed to 



have harmonic time dependence, 

Vhom(x,t) = -iu}^(x)e~ tult , 



(25) 



where uj is the oscillation frequency and £(x) is a 
function that absorbs the spatial dependence of 



Whom- Inserting this ansatz into Equation (24) 
yields the usual eigenvalue problem, 



<C4(x) = -«?£(x), 



(26) 



which, subject to appropriate boundary condi- 
tions, can be numerically solved by stellar pul- 



sation codes (e.g., Christensen-Dalsgaard 2008) 



The eigenfunctions {£jv(a;)} of C can be shown to 
be orthogonal with respect to the inner product 
defined by 



(X,Y) 



[ d 3 xp X*(x)-Y(x), 
Jv 



(27) 



see, e.g., Unno et al. (1989). Here * denotes com- 
plex conjugation and V denotes the volume of the 
star under consideration. Furthermore, the differ- 
ential operator C is Hermitian ( Unno et al.||1989" ) 
and it is generally assumed that its eigenfunctions 



form a complete set (cf., e.g., Aizenman & Smey- 
1977). As shown in Paper I, toroidal eigen- 



ers 



modes do not couple to external GWs and ad- 
ditionally these modes are of no relevance to as- 
teroseismology, since they do not perturb the hy- 
drostatic equilibrium and therefore have eigenfre- 
quencies w = ( Aizenman &: Smeyers|1977 ) . This 
allows us to solely focus on the usual spheroidal 
eigenmodes, 

£nlm( r i®><i>) = Kr,n*( r ) e r + £h,nl (r)rV\ Y lm (6, </>), 

(28) 

which we normalize according to 



JNN' 



(29) 



where N = {nlm} is an abridged index and 7 is a 
constant with dimension of a mass. We arbitrarily 
normalize the eigenfunctions to unit mass in cgs 
units, that is, I = 1 g. 

The previous discussion shows that the velocity 
field in Equation (p3| can be expanded as 



v(x,t) = ^2(-iuj N )A N (t)£ N (x)e~~ 



(30) 
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where the time-dependent complex expansion co- 
efficients A]\r(t) are called amplitudes and where 



only eigenfunctions of the form ( 28 ) need be con- 
sidered. Furthermore, we assume weak damping, 
i.e., r/N <C ujni where 77 at is the damping rate. 
Thus, we have \dAjy/dt\ <C wjv|Ajv|. In order to 
obtain analytical expressions for the amplitudes 
An, we proceed as in Paper I and solve Equation 



(231 for each velocity eigenmode, 



v N (x,t) = -iu N A N (t)£ N (x)e 



—iujj^t 



(31) 



separately. Inserting the ansatz (31 ) into Equation 



( 23 1 and adopting the weak-damping assumption 



yields the expression 



Po£ 



dA 



N 



N~ 



dt 



V' N (£ N )A N = ' 



df 



GW 



2ujn 



2w% 



dt 



The differential damping operator T> may in gen- 
eral contain time derivatives (cf. Equation (13) 



of Samadi & Goupil 2001 ) and thus the resulting 



operator V' N after applying the weak-damping as- 
sumption and factoring out ANe~ lu)Nt depends on 
the particular eigenmode N through the eigenfre- 



quency ujn- Multiplying Equation (32 1 by £* N , in 



tegrating over the volume V of the star, employing 



the normalization condition (29), and rearranging 
yields 



dA N _ e ivNt 
—j—+cr N A N = 
at 



2uj 2 n I 



d 3 x£ 



Of 



GW 



dt 



, (33) 



where we defined 



d 3 xCN ,,D ' N (£ N ) = Vd,n + i5ui 



N • 



(34) 

It is worthwhile to note that the damping oper- 
ator T> only includes dynamical sources of damp- 
ing, since the mode energy equation is not con- 
sidered. Therefore, the real part t]d,n of <7n in 
Equation ( 33 1 will contribute to dynamical linear 



mode damping. This is why we will henceforth re- 
place T)d,n by a total damping rate tjn to include 
possible additional damping mechanisms not yet 
considered in Equation ( 23 ) . The damping op- 



erator gives also rise to a dynamical shift of the 
oscillation frequency, which is manifested in the 
imaginary part Slon of on- These frequency shifts 
are generally expected to be smaller or on the or- 
der of the frequency shifts that would result from 



a non-adiabatic calculation (cf. Samadi & Goupil 
2001[ ) and hence they need not be considered here. 
Consequently, we set on — Vn- 



Equation ( 33 ) is of the same form as the corre- 



sponding Equation (23) in Paper I. Thus we can 
apply the same procedure as presented in Section 
5 and Appendix A of Paper I to solve this dif- 
ferential equation, thereby obtaining the following 
general result for the mean-square amplitude: 



(\A N (t)\ 2 ) = 



dt e 2 "« (*"-*) 



R 8 xl 

16u 2 N P 

c2(t-t ) 

I dr (K%(r))e- 

-2(t-* ) 



(35) 



where R denotes the radius of the star and () 
denotes an ensemble average, i.e., an expecta- 
tion value. In solving Equation ( 33 1 explicit use 



of the weak-damping assumption and the long- 
wavelength approximation is made. Furthermore, 
it is assumed that the driving term / GW is "turned 
on" at t = — 00 and that one is interested in the 
case of large t, i.e., in the late-time behavior of the 
system, rather than in its transient behavior. It is 
important to note that, as shown in Appendix A of 
Paper I, An is only non- vanishing for quadrupo- 
lar (I = 2) eigenmodes, that is, in Equation (35) 
we have N = {n2m} and {\A N {t)\ 2 ) = other- 
wise. Furthermore, in Equation (35) some of the 
definitions of Paper I have been transferred and 
adapted to the present hydrodynamic formalism: 
the quantity x„ is defined by the following inte- 
gral: 

Xn= / dzp (z)z 3 [£ r ,n2(z) + 3^, n2 (^)] , (36) 

Jo 

and moreover, we defined the correlation function 
K% (r) = h m (to -I) Kn (to + I) , (37) 
where h m denote the spherical components of the 



GW tensor (a definition is given by Equation ( B2 ) 
in Appendix |b|) . 

If the surrounding GW field is stationary, that 
is, all statistical quantities of the GW field are time 
independent and solely depend upon time differ- 
ences, the amplitudes become time independent, 
thanks to (r) being independent of to in this 
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case. As shown in Paper I, the mean-square am- 
plitudes are then given by the following compara- 
tively simple expression: 



(l^vl 2 ) 



R s xl 



Z2n N uj 2 N I 2 



(38) 



The second factor in this formula, P m {uS), is the 
power spectral density of the spherical accelera- 
tions h m of the GW tensor, 

/oo 
dr (K m {T)) e- iuT , (39) 
-oo 

and thus this factor solely depends on the sur- 
rounding GW field. It is worth noting that, by 
contrast, the first factor on the right-hand side of 



Equation (38) only depends on the stellar model 
and quantities derived therefrom. Consequently, 
the latter quantity describes the star's general 
ability to resonate in a particular eigenmode N 
when the star is excited by GWs, and this is why 
we christen it resonance factor. 

As in Paper I, we define the corresponding in- 
trinsic mean-square surface velocity of a particular 
eigenmode by 



,J N,s 



dflv N (x,t)-v* N (x,t) )(R) 



y N (\A N (t)\ 2 )*%(R), 



where 



(40) 



(41) 



with 1(1 + 1) = 6, () denotes time average and f2 
denotes the solid angle. We note that the defini- 
tion ( |40| docs not take visibility factors into ac- 
count, such as geometrical and limb darkening ef- 
fects, and it is equivalent to the one employed by 



Bclkacem et al. (2009). Inserting Equation (38) 



into Equation ( 40 ) yields the following expression 



for the rms surface velocities: 



R 8 xl 

6in N I 2 



P m (L> N )V 2 N (R) 



(42) 



Finally, it is interesting to note that the square 
amplitudes as well as the mean-square surface ve- 
locities show a characteristic independence on 
the radius of the star. However, this sensitive de- 
pendence on the radius is reduced to an effective 
i? 2 -dependence due to the adopted normalization 



scheme (cf. Equation (29)) 



3. Sources of gravitational radiation 

In this section, the power spectral density 



P m (w) defined in Equation (39) is calculated for 
two general types of GW radiation: radiation from 
astrophysical sources that can be individually "re- 
solved" and that emit approximately monochro- 
matic GW radiation, and an SBGW of astrophys- 
ical or cosmological origin. We start with dis- 
cussing the former case. 

3.1. Excitation by distinct astrophysical 
sources 

In the following, it is assumed that the partic- 
ular source that emits approximately monochro- 
matic radiation is sufficiently far away from the 
star under consideration, such that the incoming 
gravitational radiation can be appropriately de- 
scribed by a plane wave field as defined below. 
We assume that the incoming plane wave of fre- 
quency uj has its wave vector k aligned with the 
z-axis, k = \k\e z . In this case, the GW tensor in 
TT gauge (cf. Equation Q) takes the form 



hij(t) 




(43) 



where 
Section 



the long-wavelength approximation (cf. 
2.1 ) has already been adopted. The polar- 
e 2 ) are defined in 



ization tensors e?Je z ) and e 



Appendix [B] In order to make this formalism di- 
rectly compatible with waveforms of gravitational 
radiation from binary systems, we assume that 
A+ and A x are in general complex numbers. This 
is an important detail, which we make use of in 
Section |4.5| Finally, we adopt the usual conven- 
tion for the GW tensor that the real part is taken 
at the end of the calculation. 

The spherical components h m of hij , defined by 



h m yf 



(44) 



(cf. Append ix [b] ), are determined with the help 
of Equation ( B9 ) , where is given by Equation 
(43). One immediately finds: ho = 0, h±% = 0, 



and 



'±2 



-{-uj 2 )(A + TiA x )e- 



(45) 
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Consequently, gravitational radiation propagating 
along the z-direction solely couples to internal nor- 
mal modes of the star having an azimuthal order 
of m = ±2. 

Knowledge of h±2 now allows for the compu- 
tation of the correlation function K^It). Insert- 



ing Equation ( 45 ) into Equation ( 37 1 we find after 



some elementary transformations 
8tt 



15 



7 UJ 



± 2Re(iA + A* 



UT . (46) 

We note the high sensitivity of |-fC±2( T )l *° the 
frequency u> given by a a; 4 power law. Finally, by 
Fourier-transforming if^_° 2 (r) in accordance with 
Equation ( 39 ), one obtains the power spectral den- 
sity of h m , 



P±2{un) 



16tt 2 



-co 



I A, 



±2Re{iA+A* x )]S(oj -oj n ). (47) 



As can be read off Equation ( 46 ) , the GW field 
is indeed stationary, i.e., the correlation function 



is independent of to, and thus Equations (38) and 



( 42 1 apply. Resulting square amplitudes and in- 



trinsic rms surface velocities of excited normal 
modes (N — {n, 2, ±2}) are, respectively, given 

by 



1-4 



N 



2 E>8 2 2 



1-4 > 



±2Re(iA + A* x )] (48) 



and 



R 8 x 



2, ,4 



N 



l-4> 



(49) 



60 r] N P 

±2Re(iA + A* x )]V 2 N (R), 

where for each mode N it was assumed that there 
is a GW source at the same frequency, lo = ojn- 



Numerical results for Equations (48) and (49) in 



the case of the Sun and in view of present ampli- 
tude estimates for known and expected individual 
astrophysical GW sources are presented and dis- 
cussed in Section |4~B1 



3.2. Excitation by a stochastic background 
of gravitational waves 

In this subsection, an arbitrary SBGW is con- 
sidered and no assumption about its origin is made 



at this point, that is, we do not yet focus on a spe- 
cific background, such as a background due to cos- 
mic strings or due to astrophysical sources. This 
step is postponed to Section [4.6. 1| 

It is convenient to characterize a stochastic 
background by the dimensionless function 



0, GW (u) 



1 dp GW 
p r d In v 



(50) 



which measures the energy density of GWs per 
unit logarithmic interval of frequency in units of 
the present critical energy density, p c , that is 
needed to close the universe, 



3c 2 Hi 
8ttG ' 



(51) 



Here Hq denotes the present Hubble expansion 
rate, c the speed of light, and G the gravitational 
constant. In the following, we assume that the 
background is 

(i) isotropic, 

(ii) unpolarized, 

(iii) stationary, and 

(iv) Gaussian. 

Each of these assumptions is well motivated in the 



case of an SBGW (Allen & Romano 1999 Mag- 



giore 2008). In particular, the stationarity prop- 



erty of the background is important, since it al- 
lows us to apply the formulae for the mean-square 
amplitudes and rms surface velocities derived in 



Section^ Equations (|38j) and (|42j). Given the 
assumptions (i)-(iv), one can show that the prop- 
erties of an SBGW can be completely specified by 



a single function Sh{v) (Allen & Romano 1999 
Maggiore||2008[ ), defined by 



{h* A {v,n)h,A'{v' ,n')) 

^Siv-v'^in-n'^, 8 ^, (52) 

8-7T 

where S 2 (n — n') = S((f> — </>')6(cos6 — cosO'). 
Moreover, Iia{v) are the Fourier amplitudes of the 



GW tensor hij defined in Equation (B17), with 



A = {+, x} labeling the two polarization states, 
and () again denotes an ensemble average (i.e., 
the expectation value). The spectral density of 
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the stochastic background Sh{v) is a real, non- 
negative function satisfying Shiy) — Sh(~v) and 
it can be shown that the following relation holds: 



Ogw(^) 



in 2 
3Hl 



(53) 



(cf., e.g., |Allen fc Romano||1999||Maggiore||2008 |. 

In order to calculate the power spectral density 
P m (uS) for an SBGW, we start with the correlation 
function K^(t), 



:k '» : 1 r): " - ^hWktf)) (yfrY ytr 



dvdv' / dVtdW 



225 A A, 

A, A' 

x(/ lj4 ( l /,n)/ ) 4,(^,n'))e^(n)e^'(n') 

x^v^-^v 2 ™' 4 ' {y 2 p)*ylr, 

(54) 



where in the first step Equations (37) and ( |B9[ ) 
were employed. Moreover, we set t = to — r/2 and 
t! = to + t/2. In the second step, the plane wave 
expansion for the GW tensor of the SBGW in the 
long-wavelength approximation was used, 

M*)=E / dv [ dnh A (u,n)ef j (n)e- i27rvt , 

A J — CO J 

(55) 

cf. Equation (B18l. Inserting Equation (52) into 



Equation ( 54 ) and eliminating the delta functions 



by integrating over v' and f2' yields 

{k%{t)) = ^ J jvv"s h {vy 2 ^ J dn 
*EEE e »4(«) {yfpyy 

A i,j k,l 



2m 
kl ■ 



(56) 

Clearly, the right-hand side does not depend on to 
and thus the correlation function solely depends 
upon time differences. Laborious calculations us- 



ing Equations (|B3j)-(|B7j) and (|BL2|h(|B15] show 
that the second factor on the right-hand side of 



Equation (56) 



^EEE e » e «w WPYytr, 

A i,j k,l 

(57) 



is a constant that is independent of the azimuthal 
order m, 



6 Vm = {0,±1,±2}. 



(58) 



It is interesting to note that for m = cross- 
polarization components do not contribute to 
the expectation value of the correlation function. 
However, in the case of in = {±1,±2} the con- 
tributions of both polarization components are 
roughly of the same order. 

Consequently, one is left with 



(K£(r)> = 



75 



dwu A S h (u)e iUT , (59) 



which after Fourier-transforming in accordance 



with Equation (|39| yields the desired expression, 

(60) 



75 



or, equivalently, expressed in terms of Ogw(w) us- 
ing relation ( 53 ) , 



P m {f>) 



32ti- 

~2fT 



(61) 



Equations (60) and (61 1 nicely reflect the assump- 
tions (i)-(iv) and it should be noted that P m (w) 
is independent of the Hubble expansion rate Ho, 
so that experimental uncertainties related to Ho 
do not enter the present formalism (cf. Equations 
(|50j) and ([51])). 

However, in some cases it might be more conve- 
nient to have an expression available that directly 
relates P m (u)) to the characteristic strain spectrum 
of an SBGW. Under the assumptions (i)-(iii), the 
characteristic strain spectrum is related to Oqw 



by the following equation ( 


Allen 


1997 


2000 


Depies & Hogan||2007 


1: 



^gwM = -^u 2 h c {u) 2 . 

Consequently, when substituted we have 
16tt 2 



P m {u) 



75 



-oj 6 h c (wY 



(62) 



(63) 



Up to this point, the theory developed so far is 
still very general. In order to obtain numerical 
results for the mean-square amplitudes and rms 
surface velocities of normal modes excited by an 
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SBGW, one now has to either specify a particular 
spectral energy density r2cw(w) or, equivalently, 
a particular characteristic strain spectrum h c (u>) 2 , 
thereby focusing on a particular stochastic back- 



ground. This is done in Section 4.6.1 



Finally, using Equation (61) we reformulate 



Equations (38) and H21 in terms of the spectral 



energy density Qgw(l}) to obtain the general ex- 
pressions in the case of an SBGW, 



(l^vl 2 ) 



25 t]niunI 2 



(64) 



and 



VN,. 



R Xn WAT 



50 I 2 



Vn 



H 2 n GW (uj N )y 2 N (R) 



(65) 

As in Section|3.1[ the numerical evaluation of these 



formulae for the mean-square amplitudes and rms 
surface velocities is postponed to the next section. 

4. Numerical results for the Sun 

The theory developed in the preceding sec- 
tions is very general and applies to any non- 
relativistic star and to arbitrary SBGWs or (quasi- 
monochromatic GW sources. In this section, we 
focus on the Sun as our nearest (non-relativistic) 
star and present numerical results for resonance 
factors, mean-square amplitudes and rms surface 
velocities by making additional assumptions on 
the incoming gravitational radiation. We start 
with discussing the quantities that are directly re- 
lated to or derived from the chosen solar model. 

4.1. Eigenfrequencies and eigenfunctions 

In order to compute resonance factors, mean- 
square amplitudes, and surface velocities, knowl- 
edge of the eigenfrequencies ujn and the radial 
as well as horizontal eigenfunctions, £, r ,n2(f) and 
£/i,ra2(?"), °f I — 2 solar g- and p-mode oscillations 
is required. We note again that in the frame- 
work of the linearized theory of general relativity 
solely quadrupolar eigenmodes of non-relativistic 
stars can be excited by external GWs (cf. Sec- 
tion 2.3). Eigenfrequencies and eigenfunctions are 



computed using the standard Aarhus adiabatic os- 



cillation package ADIPLS ( Christensen-Dalsgaard 
|2008[ ) and the underlying solar model that we em- 
ploy is Model S of Christensen-Dalsgaard et al. 



-5 
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100 



1000 



10000 



v [(LiHz 



Fig. 1. — Combined damping rates of I = 2 solar g 
and p modes (from left to right) . Damping rates of 



g modes are taken from |Belkacem et al. (2009); p 



mode damping rates are extracted from MDI data 
(see the text for details). 



(1996), which uses the OPAL opacities and the 
OPAL equation of state. Eigenfunctions are nor- 



malized according to Equation (29). 



4.2. Damping rates 

As can be seen from Equation ([38] ) , or, alterna- 
tively, from Equations (48) and (64), square am- 



plitudes are the result of a balance between driv- 
ing by external GWs and damping processes in 
the interior of the Sun and thus, damping rates 
are important quantities to be reliably specified. 
In the case of solar p modes, we employ current 
experimental values obtained from data that has 
been recorded with the Michelson Doppler Imager 
(MDI) on board the SOHO spacecraft (J. Schou, 
2006, private communication). The crucial point, 
however, is the computation of g-mode damping 
rates, since they cannot be obtained from obser- 
vational data. 

Damping rates for solar g modes are taken from 
Belkacem et al. (2009). These values have been 



computed with the non-adiabatic pulsation code 



MAD (Dupret 20021, in which a time-dependent 



convection treatment presented and discussed in 
Grigahcene et al. (2005) is implemented. In Fig- 
ure^ the combined damping rates of solar p and g 
modes are depicted. Both data sets have been in- 
terpolated onto the adiabatic frequency grid com- 
puted with the ADIPLS code. 
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A few comments on (/-mode clamping rates are 
in order. For v < 110 /iHz, 77 is a decreasing func- 
tion of frequency and follows a very regular law, 
77 cx v^ 3 , which can also be recovered from quasi- 



adiabatic analytical computations (sec Belkacem 
et al.||200"9 ). This result is due to radiative damp- 
ing in the inner core being the dominant contri- 
bution to the work integral and thus to the damp- 
ing rates (Belkacem et al. 2009). However, for 
77 > 110 /iHz, time-dependent convection terms 
become significant and damping rates are sensi- 
tive to the parameter (3 (Grigahcene et al. 2005) 



in the perturbation of the energy closure equa- 
tion. This parameter is used by Belkacem et al. 
( 2009 ) to model convection-pulsation interactions 



and the value adopted by these authors was cho- 
sen such that good agreement between theoreti- 
cal and observed damping rates of solar p modes 
is achieved. However, there is no evidence that 
this value can readily be used in (/-mode calcula- 
tions. Consequently, Belkacem et al. (20091 con- 



clude that (/-mode damping rates for 77 > 110 /iHz 
are questionable. However, uncertainties related 
to the modeling of convection-oscillation interac- 
tion are not inherent in (/-mode damping rates be- 
low 110 /iHz. 

4.3. Long-wavelength approximation 

The long-wavelength approximation (see Sec- 
if appropriate will apply to all GW 



2.1) 



tion 

sources, since it only requires the quotient Igw = 
gJ/(A/2tt) = (2irdiy)/c to be very small compared 
to unity. Here, d is the diameter of the star un- 
der consideration, c is the speed of light, A de- 
notes the wavelength of the incident GW radia- 
tion, and v its frequency. In the solar case one 
has Zqw < 6 x 10~ 3 for the (/-mode frequencies 
(y < 200 /tT1z) and Zgw < 0.15 for the p-mode fre- 
quencies (y < 5 mHz) considered here and thus the 
long-wavelength approximation is well justified in 
the solar case. For stars more massive than the 
Sun radii increase whereas (/-mode oscillation fre- 
quencies decrease. Generally, the long-wavelength 
approximation is expected to be appropriate for g- 
mode frequencies of all non-relativistic stars and 
it is also expected to be well motivated for p-mode 
frequencies. However, the smallness of Zqw has to 
be checked case by case. 




Fig. 2. — Numerical values for Xn (top) and for 
resonance factors (bottom) of / = 2 solar g and p 
modes (from left to right), g modes show a much 
higher general ability for resonance in the case of 
external forcing by gravitational waves. 

4.4. Resonance factors 



In Section 12. 31 we introduced the resonance fac- 



tor, 



Qr 



R Xn 

cs.N = ™ FTTo- (66) 

It contains all quantities in the expression for the 



square amplitudes, Equation (38), that are solely 



dependent on the solar model and that do not de- 
pend on the incoming GW radiation. Hence, the 
resonance factor associated with a particular nor- 
mal mode provides a general measure for the abil- 
ity of the Sun to resonate in this particular eigen- 
mode, independent of the assumed external GW 
field. 

Numerical results for the resonance factors 
Qres,N based on the computations described in 
Sections |4.1| and |4.2| are depicted in Figure [2] 
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In addition, we present numerical values for the 



quantity Xn i which is defined in Equation ( 36 ) 



and which is among the main contributors to the 
resonance factor. We note that the high values 
for the resonant factors present in Figure [2] result 
from the normalization of the eigenf unctions (cf. 
Section 2.3). Only relative values have a direct 



meaning and are to be considered. 

From Figure [2j it is evident that the overall 
spectral shape of Q res ,N is dominantly determined 
by Xn > with the offset between the g-mode and the 
p-mode branch being due to division by damping 
rates and squared frequency. While a power law 
would only be a crude first-order estimate for p 
modes, it is interesting to note that the g-mode 
resonance factors very regularly follow a power 
law. The latter fact is due to the power-law behav- 
ior of the damping rates in the g-mode frequency 
range (cf. Section 4.2), together with the power 
law for g-mode frequencies present in the numeri- 
cal results for the quantity Xn (top plot of Figure 

As is also obvious from Figure [2] g modes show 
up to seven orders of magnitude higher resonance 
factors than p modes and thus g modes have a 
much higher general ability to resonate in the case 
of this excitation mechanism. However, in order to 
be able to decide whether or not solar g modes are 
indeed excited more intensively than p modes, one 
has to compare the corresponding mean-square 
amplitudes. 
and l4lQl 

4.5. Excitation by distinct astrophysical 
sources 

In this section, we present numerical results for 
the mean-square amplitude and rms surface ve- 
locity of a particular solar normal mode that is 
excited by a quasi-monochromatic and resolvable 
astrophysical source, by which we mean that the 
number of sources in the same, sufficiently wide 
frequency bin does not exceed one, or, at least, 
the GW signals of other sources in this bin should 
be comparatively small, such that they can be 
neglected. Apart from an SBGW of cosmologi- 



This is addressed in Sections 4.5.2 



cal origin (see Section 4.6), ultra-compact binary 
systems are among the most relevant astrophysi- 
cal sources for gravitational radiation in the low- 
frequency band of stellar oscillations to be "ob- 
served" from the solar system. In particular, AM 



CVn stars, which are interacting ultra-compact 
binaries (white dwarfs accreting helium-rich ma- 
terial from a compact companion), are the only 
known sources in this frequency range (see Section 
4.5.11, albeit many other astrophysical sources 



like massive black hole mergers, compact objects 
spiraling into massive black holes, etc., are ex- 
pected. In view of this it seems advisable to ex- 
plicitly evaluate the still very general formalism 
presented in Section |3.1| for the case of a binary 
system. In this section, we are only interested 
in the quasi-monochromatic radiation associated 
with binaries before the final spiral-in phase, that 
is, the timescale for frequency shift due to loss 
of energy by gravitational radiation is large com- 
pared to the GW period. We review the expected 
intrinsic strain amplitudes for known and expected 
systems. This allows us to deduce a realistic order 
of magnitude estimate for the expected strain am- 
plitudes to be employed in calculations for square 
amplitudes and surface velocities; the results of 



these computations are presented in Section 4.5.2 



4-5.1. Assumptions on the incoming gravita- 
tional radiation 

According to the quadrupole formula, gravita- 
tional radiation in TT gauge from a binary system 



in circular orbit is given by ( Timpano et al.||2006 
Maggiore|[2008l ) 



h^t) = A+e-^e+fe^+Axe-^e^e,), 

(67) 

where the polarization amplitudes are 
4 



, Lj = _ (GM ) 5 / 3 (™ gw ) 2 / 3 
1 + cos 2 



re* 



,i27ri/ gw r/c 



(68) 



A_=i-^(GM) 5 / 3 (nv sw ) 2 / 3 

xcosee t2 ™^ /c . (69) 

Here M = {M x M 2 f /h /(Mi + Af 2 ) 1/5 is the chirp 
mass, with M\ and M 2 being the component 
masses, O is the inclination of the binary, r its 
distance, and z/ gw denotes twice the instantaneous 
orbital frequency. Moreover, we already assumed 
that the radiation is coming from the —e z direc- 
tion and we also adopted the convention intro- 
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duced in Section 3.1 of allowing complex quan- 
tities and taking the real part at the end of the 
calculation. We note that adopting this conven- 
tion allows us to make the waveforms of a binary 
system directly compatible to the form of gravita- 
tional radiation assumed in Section 13.11 More- 
over, we note that galactic double white dwarf 
(DWD) systems, DWD systems being the most 
common binary systems, are expected to have in- 



mans 


2009 


Ruiter et al. 


2010) 



In addition, the 
assumption of a circular orbit is also very well sat- 
isfied in the case of galactic AM CVn stars, which 



are also commonly expected to be circular (Nele- 


mans et al.|2004 Roelofs et al.|2006||Roelofs et al. 


2007a 


Ruiter et al.||2010). 


Following 


Timpano et al. (2006 


), Roelofs et al. 


(2006) 


, and 


Roelofs et al. (2007a 


, we define the 



GW strain amplitude h by 



h = 



\A- 



1/2 



(70) 



Furthermore, we take the mean power spectral 
densities P±2(ujn) of N = {rt,2,±2} modes as 
a measure for predicting the effect of such indi- 
vidual astrophysical sources on solar oscillations, 
such that the corresponding mean-square ampli- 
tudes and rms surface velocities (cf. Equations 



(481 and (49)) become 



1-4 



and 



N 



r^2 



7H R*X>N h 2 

15 m P ' 



VN,s 



30 r, N P 



N 



(R) 



(71) 



(72) 



with the bar over P± 2 (ujn), \An\ 2 , and vjsr iS indi- 
cating that the mean of the N = {n, 2, ±2} modes 
is considered. 

The only known, guaranteed GW sources in 
the low-frequency band that have well-determined 
strain amplitudes with reliable error bars are five 
AM CVn systems, which are listed in Table [T] 
(Roelofs et al. 2007a). These sources have strain 
amplitudes of the order of h ~ 10~ 22 . In addi- 
tion, there are several tens of known DWD bina- 
7 x 10~ 7 



ries (v. 



4 x 10- 4 Hz; cf. Nele- 



mans et al. 2005) and AM CVn systems (u g 



Table 1: Predictied Gravitational Wave Strain 
Amplitudes h and Frequencies ^ gw of Five AM 



CV n Binary Systems (|Roelofs et aL]|2007a[) 



Star 




h 






(mHz) 






AM CVn 


1.944 


2-0 +U o1 x 


10 -22 


HP Lib 


1.813 


o 7 +0.6 


1Q -22 


CR Boo 


1.360 


2 l +0 - 4 x 

z - x -0.5 x 


1Q -22 


V803 Cen 


1.253 


3.0±°- 5 7 x 


10 -22 


GP Com 


0.7158 


[4.0 - 6.6] 


X 10~ 23 



6 x 10 3 Hzj 3 } t heir strain amplitudes 



Nelemans et al. 2004), however, can 



5 x 10~ 4 - 
(cf., e.g., 

only be estimated for lack of accurate distance 
measurements. Moreover, focusing on frequencies 
v gw ~ 10~ 3 — 10~ 2 Hz, more than 10,000 galactic 
DWD binaries and thousand to more than 10,000 
AM CVn systems are expected in the strain ampli- 
tude range of h ~ 10~ 24 — 10~ 22 (Nelemans et al 
2004} [Roelofs et al.||2007b[ |Ruiter et al.||2010[ ). 

Besides these known sources, massive black hole 
binaries in the mass range 10 3 — 10 s Mq are ex- 
pected to provide even stronger gravitational radi- 
ation at /iHz and mHz frequencies. For example, 
a representative circular massive black hole binary 
with masses M 2 = O.LMi = 10 5 M Q at redshift 
z = 1 produces a strain amplitude of h ~ 10 -20 
at i/gw ~ 10 -5 Hz before the final spiral-in phase 



( |Sesana et aT1|2005| ). 

The preceding discussion shows that for an or- 
der of magnitude estimate for amplitudes and 
surface velocities of solar normal modes excited 
by a particular quasi-monochromatic astrophysi- 
cal source, it is reasonable to assume the source to 
have maximal strain amplitudes of h ~ 10 -20 . 

4-5.2. Numerical Results 

In order to provide an order of magnitude es- 
timate for the effect of external forcing by quasi- 
monochromatic GW sources on solar oscillations, 
we assume that for each eigenfrequency of the 
Sun there is a GW source at the very same fre- 
quency with a gravitational strain amplitude of 



3 A list of known AM CVn stars 
ties, and corresponding references 



their proper- 
available at: 

http : / /www. astro . ru.nl/$\sim$nelemans/dokuwiki/ 
\doku .php?id$=$verif ication_binaries : am_cvn_stars 



13 




100 1000 
v frtHz] 



10000 



Fig. 3. — Numerical results for mean-square am- 
plitudes (top) and rms surface velocities (bottom) 
of I = 2 solar g and p modes (from left to right), 
assuming a constant gravitational wave strain am- 
plitude of h = 1 x 10~ 20 over frequency from in- 
dividual astrophysical sources as discussed in Sec- 
tions ED and HSU 



h = 1 x 10~ 20 (cf. Section 4.5.1 ). Using the damp- 



ing rates discussed in Section 4.2 eigenfunctions 
and eigenfrequencies as well as quantities directly 
related to the solar model, the former and the lat- 
ter being discussed in Section |4~T| and Equations 
( 71 ) and |72|) , we obtain the results depicted in 



Figure[3] 

Again we note that in the plot showing the 
mean-square amplitudes, only relative values have 
a direct meaning, since the high absolute values 
are an immediate consequence of the adopted nor- 
malization condition for the eigenfunctions (cf. 
Section 2.3 1. 



As is obvious from Figure [3j low-frequency g 
modes show smaller mean-square amplitudes than 



low-frequency p modes. This is due to suppression 
by the w 4 dependence on frequency present in the 



power spectral density P m (uj) (cf. Equation (47) 
and the results for the resonance factors in Fig- 
ure]^]). However, mean-square amplitudes of high- 
frequency g modes (v gw > 80 /xHz) exceed the cor- 
responding values of low-frequency p modes by up 
to five orders of magnitude. This shows that, due 
to the definition of the square amplitudes, Equa- 



tion ( 30 1 , and due to the normalization condition 



( 29 1 , high-frequency g modes are more strongly 



excited than p modes by several orders of magni- 
tude. 

In Figure [3j we also show numerical results for 
rms surface velocities, which are found to reach 
values of vn, s ^5 10 _8 mms~ 1 . With perhaps a 
few exceptions, g-mode values generally lie be- 
low p-mode values. Consequently, although high- 
frequency g modes are excited more strongly than 
p modes, this behavior is not obvious any more 
at the surface. This is because of g modes being 
exponentially damped in the convection zone. In 
the analytical expressions for the rms surface ve- 



locities, Equations (42 1, (49), and (72), damping 



of g modes in the convection zone can be seen to be 
modeled by the quantity ^j^(R), which is defined 



in Equation ( 41 ) and which essentially measures 



the amplitude of the eigenfunctions at the surface 
of the star. Numerical results for this quantity 
are depicted in Figure |4j in which damping of g 
modes in the convection zone is nicely reflected by 
<?-mode values being orders of magnitude smaller 
than p-mode values. 

We finally note that both mean-square ampli- 
tudes and rms surface velocities for g modes show 
very regular power-law behavior. This is due to 
the power-law dependence on frequency shown by 



the damping rates (cf. Section 4.2), the quantity 



Xn (cf. Section 4.4), and the power spectral den- 



sity P m (cf. Equation (47)). In the case of surface 



velocities, it is important to note that the quan- 
tity ^^(i?) also shows a power law for asymp- 
totic (low-frequency) g modes. The deviation from 
power-law behavior inherent in rms surface veloc- 
ity values at high g-mode frequencies is predomi- 
nantly due to a departure from power-law behav- 
ior of ty%(R) at frequencies above 100 Hz, which 
is evident from Figure [2] Assessing the signifi- 
cance of the surface velocities found in this subsec- 
tion by comparing them to corresponding values 
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Fig. 4. — Numerical values for the quantity ^ 2 N (R) 
(cf. Equation (41)) of I — 2 solar g and p modes 
(from left to right); g-mode values are orders of 
magnitude smaller than p-mode values. A power- 
law behavior for asymptotic g modes is clearly vis- 
ible. 



in the case of stochastic excitation by turbulent 
convection (first two terms on the right-hand side 
of Equation |20|) is postponed to Section 4.7 



4.6. Excitation by a stochastic background 
of gravitational waves 

In this section, we discuss numerical results for 
mean-square amplitudes and rms surface velocities 
of solar normal modes excited by an SBGW. As is 



pointed out in Section 3.2 conducting such calcu- 
lations requires to either specify a particular spec- 
tral energy density Q gw (uj) or, equivalently, a char- 
acteristic strain spectrum h c (uj) 2 for the SBGW to 
be expected. Therefore, we start with briefly re- 
viewing possible sources for an SBGW in order to 
deduce a realistic assumption on the spectral en- 
ergy density to be employed. 

4-6.1. Assumptions on the incoming gravita- 
tional radiation 

A (rather weak) contribution to the SBGW 
at /LtHz frequencies is expected from astrophysi- 
cal sources. The characteristic strain spectrum 
of an SBGW originating from coalescing massive 
and supermassive black hole binaries at the cen- 
ter of galaxies is expected to peak at frequencies 
- 10 Hz, with hi 



'gw 



10" 



10 



" 28 ( Sesana et al 
2004[ |Enoki et al.||2004[ ); at v gw ~ 10" 4 Hz these 
authors find hi ~ 10 -38 which translates into 



JIqw ~ 10~ 10 using Equation (62 1. The upper 
bound on the entire frequency range of 10 -12 Hz < 
^gw ~ 10 _1 Hz is Ogw ~ 10~ 8 , which is reached at 
v gw ~ 10~ 6 Hz. A background due to galactic and 
extragalactic detached white dwarf binaries is gen- 
erally weaker than the one by black hole binaries 



and has been calculated by Nelemans et al. (2001) 



and Farmer & Phinney ( 2003 ) , respectively. In the 



case of galactic binaries, the former authors find 
JIqw ~ 10~ n at i/ gw ~ 10~ 4 Hz, with an upper 
bound of f^GW ~ 10~ 10 on the whole frequency 
range, 10 -4 Hz < z/ gw < 10 -3 Hz, whereas the lat- 
ter authors in the case of extragalactic binaries 
obtain ft GW - 10~ 13 - 10~ 12 at v gw ~ 10~ 4 Hz, 
with an upper limit of flaw ~ 10~ n on the entire 
frequency range, 10 -6 Hz < i/ gw < 10 _1 Hz. 

The second general category of an SBGW is 
formed by cosmological sources. A weak contribu- 
tion at /iHz frequencies is expected from an SBGW 
that originates from quantum fluctuations dur- 



1990 


Bar-Kana 


1994 





10" 



After a decrease with frequency as v gVl 
from 10~ 18 Hz to 10~ 16 Hz, typical inflation mod- 
els predict almost scale-invariant spectra of GWs 
over many orders of magnitude in frequency, e.g., 
10~ 16 Hz < f gw < 10 _1 Hz, up to a region near a 
cutoff frequency, which is related to the end of re- 
heating. The constant spectral energy density typ- 
ically found in this quasi-scale-invariant frequency 
regime is Oqw ;$ 10 -14 . An SBGW arising from 
amplification of vacuum fluctuations in string the- 
ory inspired pre-big-bang models is thought to be 
very weak in the low-frequency range of stellar 
<?-mode pulsations but it would be among candi- 
dates for a dominant SBGW at high frequencies, 
i/ gw > 1 Hz. Such models generally yield a spec- 
trum that shows an increase at low frequencies 
proportional to is 3 up to a turn-over frequency v s , 
above which it is found that f^cw ~ ^ 3 ~ 2ai , and 
finally the spectrum cuts off at a frequency v\, 



which is expected to be around 10 10 Hz (Brustein 


et al.||1995 


Buonanno et al.|1997 Mandic & Buo- 


nanno 


20061. Depending on the model parame- 



ters, v s , n, and fi, these models are able to pro- 
duce spectral energy densities at high frequencies, 
e.g., v gw > 10 Hz, that lie close to present observa- 
tional bounds on an SBGW. However, at frequen- 
cies of stellar g-mode pulsations below 10~ 3 Hz 
the spectral energy density Qqw is thought to be 
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smaller than the SBGW predicted by the inflation 
models discussed above by orders of magnitude. 

Cosmic strings might be the dominant source 
for an SBGW in the low-frequency range of solar 
and stellar g-mode pulsations and they presently 
seem to be the only possible source to produce 
an SBGW that lies close to current observational 
bounds at low frequencies as well as at high fre- 
quencies. Cosmic strings are one-dimensional 
topological defects that can be formed during sym- 
metry breaking phase transitions in the early uni- 
verse ( |Kibble 19761. In addition, string-theory- 



inspired inflation scenarios also produce cosmic 
strings, which are referred to as cosmic super- 



string 


s (|Jones et al. |2002| Sarangi & Tye 


2002 


Dvali & Vilenkin 2004; Jones et al.|2003| Copeland 


et al. 


2004 


). Whereas in the former case of field 



theoretic strings the reconnection probability is 
unity, p = 1, cosmic superstrings have a recon- 
nection probability less than unity, with expected 
values of 10~ 3 < p < 1 ( Jackson et al. |2005| ). 
A network of cosmic (super)strings gives rise to 
an SBGW via incoherent superposition of bursts 
from cusps and kinks on these strings. This mech- 
anism of GW production has been studied inten- 
sively in the past (e.g., |Vilenkin||1981| | Caldwell 



& Allen 


1992 


Caldwell et al. 


1996; Damour & 


Vilenkin 


2001 


2005 Siemens et al. 


2007 Olmez 


et al. 2010). After a decrease or a steep rise at 



the low-frequency end of the spectrum, depending 
on the assumed mechanism that determines the 
string length, gravitational backreaction ("small 
loops") or large-scale dynamics of the network 
("large loops"), the spectral energy density Uqw 
typically shows a scale-invariant behavior over 
many orders of ma gnitude. Approxima te analyti- 
cal calculations by 



Olmez et al. 



(2010) show that 



this constant is solely determined by two model 
parameters, the reconnection probability p and 
the dimcnsionless constant G/i that characterizes 
the gravitational interaction of strings with string 
tension /i; G denotes Newton's gravitational con- 
stant. The current observational bound on the 
string tension from cosmic microwave background 
(CMB) observations is G/j, < 6.1 x 10~ 7 {Po gosian 
et al.||2004| |Jeong & Smoot||2005[ ). In the small 
loop case, the aforementioned constant is given by 



Gfi 



(73) 



whereas in the case of large loops one finds 



flow oc 



p 



(74) 



( |6lmez et aL]|2010 |. 

In view of the spectral properties of various pos- 
sible candidates for an SBGW mentioned above 
and the observational bounds to be listed below, 
we adopt the following scenarios in order to es- 
timate mean-square amplitudes and rms surface 
velocities of solar normal modes, in particular the 
g modes, excited by an SBGW. We focus on an 
SBGW due to cosmic strings and neglect the con- 
tributions of all other possible sources, since their 
contributions to an SBGW at solar g- and p-mode 
pulsation frequencies are expected to be very small 
compared to possible spectra produced by cos- 
mic strings. Furthermore, we consider the fol- 
lowing two cases. First, we adopt an optimistic 
scenario and assume a constant spectral energy 
density over the frequency range of solar g- and 
low-frequency p-mode oscillations of f2Gw( w ) = 
1 x 10 -5 . This is approximately the case at solar 
g-mode frequencies for a small loop model of ra- 
diation from cusps with parameters p = 5 x 10 -3 , 
Gu = 10~ 7 , and e ~ 10~ 7 (cf. Siemens et al. 



2007 



Olmez et al. 



2010). We note that this 



model is perfectly compatible with all observa- 
tional bounds, such as the current pulsar limit 
at low frequencies around 10~ 8 Hz (Ijenet et al. 



2006), the direct bound around 100 Hz deduced 



from the Laser Interferometer Gravitational-wave 



Observatory (LIGO) science run S5 (Abbott et al. 
2009 ) , the bound based on CMB and matter spec- 



tra (Smith et al. 2006), the indirect bound due 
to the big bang nucleosynthesis (BBN; jCyburt 
et al. 2005), and the CMB bound on the string 
tension mentioned above. Second, we consider 
a rather pessimistic case and assume a constant 
spectral energy density over the frequency range 
of solar g- and low-frequency p-mode oscillations 
of r2cw( w ) = 1 x 10 -8 . Such a value for the 
constant part of the GW spectrum would, for 
example, result from a large loop model of su- 
perstrings with G/i = 10~ n and p = 10 _1 or 
from a large loop model of field theoretic sting s 
(p = 1) with Gfi = lQ- g (|6lmez et aT1|2010|. 



Both of these models include radiation from cusps 
and kinks, they already reach the scale-invariant 
regime at lower frequencies than the <?-mode fre- 
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quency range, and both are perfectl y consistent 
with all current observational bounds ( Olmez et al. 



2010 ). For the remainder of this paper, we employ 
H = 73kms~ 4 Mpc~ 1 for the Hubbl e par ameter, 
as is done by 



Siemens et al. 



(2007) and 



Olmc; 



et al. (2010). We note, however, that our results 
are independent of the adopted value of H and 



its experimental uncertainty (cf. Section 3.2) 



4-6.2. Numerical Results 

In Figure [5] we present numerical results for 
mean-square amplitudes and rms surface veloci- 
ties of quadrupolar solar g and p modes excited by 
an SBGW. The calculations have been performed 



using Equations ( 64 ) and ( 65 1 , eigenfunctions and 



eigenfrequencies as well as quantities directly re- 
lated to the solar model as discussed in Section 
|4.1[ damping rates as presented in Section [4~2| and 
the model scenarios for an SBGW due to cosmic 



strings described in Section 4.6.1 



As depicted in the top plot of Figure [5j mean- 
square amplitudes of g modes above 40 /iHz exceed 
corresponding values for low-frequency p modes 
by up to six orders of magnitude and thus we con- 
clude that these g modes are much more strongly 
excited than p modes (cf. the definition of square 
amplitudes, Equation (|30[ ), and the normalization 
condition, Equation (29)). The reason for the 
very low-frequency (high-radial order) g modes 
having smaller mean-square amplitudes than low- 
frequency p modes is suppression of low-frequency 
modes that results from an explicit frequency de- 
pendence of the power spectral density P m (ui), 
P m (uS) oc u (cf. Equation (61) and the results for 



the resonance factors presented in Section (4.4)) 



Note that in the case of excitation by a partic- 
ular astrophysical source, P m (iS) even includes a 
w 4 suppression factor for low-frequency modes (see 



Equation ( 47 )) and as a result of this, mean-square 
amplitudes of g modes below 80 /iHz lie below cor- 
responding values for low-frequency p modes (see 
Section |4.5.2 ) . As a rather marginal remark we 



again mention that in the top plot of Figure [5] 
only relative values have a direct meaning, since 
absolute values are strongly dependent on the ac- 



tual normalization constant I, cf. Equation (29) 



Therefore, if ficw(w) is considered to be constant, 
the actual value of Slew is irrelevant to the discus- 
sion of mean-square amplitudes and we arbitrarily 
set Hgw = 1 x 10~ 5 in this plot. 




100 1000 10000 

v [/iHz] 

Fig. 5. — Top: numerical results for mean-square 
amplitudes of I = 2 solar g and p modes (from 
left to right) assuming Ogw(^) to be constant, 
SIqw = 1 X 10~ 5 . High-frequency (low- radial or- 
der) g modes are more strongly excited than low- 
frequency p modes by up to six orders of magni- 
tude. Bottom: numerical results for rms surface 
velocities of the same modes assuming a constant 
spectral energy density of J]qw = 1 x 10 -5 (top 
curve) and ficw = 1 x 10 -8 (bottom curve). The 
dotted curves in between correspond to setting 
f^GW = 1 X 10~ 6 and Qgw = 1 x 1CT 7 , respec- 
tively. 



Again, and for the same reason as stated in 
Section |4.5.2[ mean-square amplitudes of g modes 
and rms surface velocities of asymptotic (low- 
frequency) g modes show characteristic power-law 
behavior. The deviation from a power law visible 
in the rms surface velocities of high-frequency g 
modes is predominantly due to the non-power-law 
behavior of the quantity at frequencies above 



100 Hz (cf. Figure gj and Section [45j) . 
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In Figure [5] we plot numerical results for rms 
surface velocities adopting the optimistic scenario 
(upper curve, Slew = 1 x 10~ 5 in the considered 
frequency range) and the rather pessimistic sce- 
nario (bottom curve, f^GW = 1 X 10~ 8 in the con- 
sidered frequency range), which are discussed in 
Section 14.6.11 The dotted curves in between are 
plotted to guide the eye and correspond to set- 
ting f2 GW = 1 x 10~ 6 and Hqw = 1 x 10~ 7 . As in 
the case of excitation by a particular astrophysical 



source (Section 4.5.2), we again find that most of 



the g-mode values lie below the p-mode rms sur- 
face velocities, which is due to the quantity ^f^(R) 
that in the present formalism can be considered as 
to model the exponential damping of g modes in 
the convection zone. Setting f^cw = 1 x 1CP 5 we 
obtain g-mode rms surface velocities in the range 
10 _9 rrims _1 < wjv . 



of 10 _11 mms _1 < vn. 



1 x 10 -8 values are found to be of the order 
< 10 _s mms _1 , with in- 
creasing values toward higher frequencies. In this 
context, it is important to stress that rms surface 
velocities scale with the square root of Qgw(v), 



cf. Equation (65), and thus, for example, shifting 



Hqw by two orders of magnitude only changes 
surface velocities by one order of magnitude. In 
the following section, we will focus on a compari- 
son with present estimates for g-mode rms surface 
velocities in the case of stochastic excitation by 
turbulent convection, thereby assessing the signif- 
icance of the results presented here. 

4.7. Comparison with stochastic excita- 
tion by turbulent convection 

Theoretical quantitative estimates for rms sur- 
face velocities of low-degree (I — 1,2,3) g modes 
that are stochastically excited by turbulent con- 
vection differ from each other by orders of magni- 
tude. The current range for quadrupolar (/-mode 
rms surface velocities is 10~ 3 mms -1 < vn, s < 
lmms -1 (for a review on the topic see, e.g., Ap 



pourchaux et al. 2010). These order of mag- 



nitude differences are predominantly due to the 
choice of the assumed turbulent eddy-time corre- 
lation function. Using a Lorentzian profile, [Belka-| 
cem et al. (2009) find rms surface velocities of 
the order of hjv, s ~ 0.1 — 1 mms -1 for quadrupo- 
lar g modes, whereas employing a Gaussian pro- 
file theseauthors essentially reproduce the results 



O.lmms" 1 ), who also used a Gaussian function. 
A rather crude estimate of g-mode amplitudes us- 
ing the equipartition ansatz of Goldrei ch fc Kee^| 
ley (1977) and without referring to an explicit 



mode excitation and damping model was reported 



by Gough (1985), who obtained rms surface ve- 



locities for quadrupolar g modes of the order of 
Vn,s ~ 0.1 mms -1 . 

Comparing these findings with our theoretical 
estimates for rms surface velocities of quadrupo- 
lar g modes excited by GWs of a particular astro- 



physical source of vn, 
tion 



< 10 8 mms (see Sec- 



4.5.2 ), we conclude that this excitation mech- 



anism seems not to be relevant to the excitation 
of quadrupolar g modes. However, this picture 
changes dramatically if we consider excitation by 
an SBGW. In this case and depending on the 
adopted scenario (optimistic: Qgw = 1 x 10~ 5 in 
the g-mode frequency range; pessimistic: Oqw = 
1 x 10~ 8 in the g-mode frequency range; see the 



of Kumar et al. (1996) (10 _3 mms 



i < 



vn,. 



< 



discussion in Section 4.6.2), we find maximal sur- 
face velocities in the range of 10~ 5 — 10~ 3 mms -1 
(cf . Figure [5| . Consequently, present models for 
an SBGW due to cosmic strings are able to excite 
solar g modes up to values that lie very close to 
(or possibly even within) the presently expected 
range for rms surface velocities of solar g modes 
in the case of excitation by turbulent convection. 
Thus, we conclude that once the experimental sen- 
sitivity will allow for the detection of individual 
quadrupolar g modes in the near future, and pro- 
vided that rms surface velocities due to excitation 
by turbulent convection in fact lie in the presently 
predicted range, it might be possible to measure 
the effect of an SBGW on quadrupolar solar g 
modes, or, at least, it might be possible to impose 
additional constraints on an SBGW. 

5. Conclusions 

In the present paper, a hydrodynamic formal- 
ism has been developed that describes excitation 
of normal modes of non-relativistic stars by exter- 
nally generated GWs. In particular, we focused 
on stationary GW signals and we considered exci- 
tation by individual astrophysical sources, such as 
binary systems, as well as excitation by an SBGW 
of astrophysical or cosmological origin. 

The equation of motion that governs the in- 
ternal motions of a star under external forcing 
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by GWs is derived from the field equations of 
general relativity by applying the linearized the- 
ory approach, which is valid by definition of non- 
relativistic stars, by adopting the long-wavelength 
approximation, which is valid for the Sun and 
which is expected to be valid for all non-relativistic 
stars (cf. Section 4.3), and by assuming non- 
relativistic internal motions of the star. This 
Euler-like equation is then transformed into an 
equation of motion solely for the oscillations of 
the star by following the procedure discussed by 
Samadi & Goupil (20011, thereby assuming adi- 



abatic oscillations that are decoupled from the 
turbulent medium and incompressible turbulence. 
However, in contrast to Samadi & Goupil (2001) 



we did not adopt the Cowling approximation and 
we showed that terms related to perturbations in 
the Newtonian gravitational potential can be re- 
tained as far as oscillations are concerned. As a 
by-product of applying this procedure, we obtain 
a generalized and unified theoretical description 
of the mode excitation processes due to turbulent 
convection and GWs, combined in a single con- 
sistent forced wave equation. Focusing solely on 
excitation by GWs, we were able to analytically 
solve the resulting inhomogeneous linear differen- 
tial equation, with the general solution given by an 
expression for the mean-square amplitudes of stel- 



lar normal modes, Equation (35 1 . Only quadrupo- 
lar (I = 2) normal modes are excited. In the case 
of stationary GW signals, we were able to fur- 
ther simplify the aforementioned expression and, 
as the main result of the first part of this paper, 
we obtain comparatively simple relations for the 
mean-square amplitudes and rms surface veloci- 
ties of stellar normal modes, Equations p8| and 
( [42] >. In Section [3j we derived analytical expres- 
sions for the power spectral density P m (u>) of the 
GW tensor's spherical accelerations h m in the case 
of a particular astrophysical source of GW radia- 
tion and in the case of an SBGW. 

Building on these theoretical results and focus- 
ing on the Sun as our nearest (non-relativistic) 
star, we computed and presented numerical results 
for mean-square amplitudes and rms surface veloc- 
ities of quadrupolar solar normal modes excited by 
the aforementioned types of GW radiation. We 
find a very regular power-law behavior for reso- 
nance factors, mean-square amplitudes, and rms 
surface velocities of g modes. As is shown by the 



numerical results for the resonance factors, solar 
g modes have a much higher general ability to be 
excited than p modes and indeed, as shown by the 
results for mean-square amplitudes, g modes with 
frequencies above 80 /iHz (excitation by a particu- 
lar astrophysical source) or 40 /iHz (excitation by 
an SBGW) are excited more intensively than p 
modes by many orders of magnitude. Suppression 
of low-frequency (high-radial order) g modes re- 
sults from an explicit power-law frequency depen- 
dence of the power spectral density P m (u) (cx w 4 
and oc co in the above mentioned cases, respec- 
tively) . 

In the case of excitation by individually re- 
solvable astrophysical sources, we obtain maxi- 
mal rms surface velocities of quadrupolar solar g 
modes of v^, s ~ 10 mms , assuming an in- 
trinsic GW strain amplitude of a binary system 
of h = 1 x 10~ 20 , which is a realistic (albeit opti- 
mistic) assumption at /iHz frequencies as argued 
in Section 4.5.1 Compared to corresponding re- 



sults for excitation by turbulent convection, this 
excitation mechanism due to individual astrophys- 
ical GW sources yields values that are four orders 
of magnitude smaller and thus it does not seem to 
be relevant to the excitation of solar g modes, pro- 
vided that the current range for quadrupolar so- 
lar g-mode rms surface velocities that results from 
various predictions by different theoretical mod- 
els is correct (at present: 10 _3 mms 



vn,. 



< 



1 mms , see the discussion in Section 4.7). How- 



ever, this picture changes dramatically when con- 



sidering an SBGW. As argued in Section 4.6.1 



cosmic (super)strings presently seem to be the 
only source of GW radiation at low frequencies 
that is likely to produce an SBGW close to cur- 
rent detection limits. Adopting an optimistic sce- 
nario of f^GW = 1 x 10 -5 and a rather pessimistic 
scenario of Oqw = 1 x 10~ 8 in the solar (/-mode 
frequency range, we obtain maximal rms surface 
velocities of 10~ 5 — 10~ 3 mms -1 . Maximal rms 
surface velocity values thus lie very close or might 
possibly even lie within the expected amplitude 
range for quadrupolar solar g modes excited by 
turbulent convection, which is currently believed 



to be responsible for g-mode excitation (e.g., Ap 
pourchaux et al.|2010 1 . We note that both models 



for the SBGW are in perfect agreement with all 
observational bounds on an SBGW to date (see 
Section 4.6.1). We also note that when shifting 
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^GW by several orders of magnitude, rms surface 
velocities are shifted only by half the amount of 
orders of magnitude, since the latter depend on 
Ogw only through the square root of flew- 

Hence, we conclude that once it will be possible 
to detect individual solar g modes in the near fu- 
ture, one will either be able to measure the effect 
of an SBGW due to cosmic strings on solar g-modc 
oscillations, thereby detecting GWs and using the 
Sun as a detector for GWs, or one will be able 
to deduce additional constraints on an SBGW at 
/xHz and mHz frequencies. We repeat that GWs 
are only able to excite quadrupolar (/ = 2) os- 
cillation modes and that turbulent convection ex- 
cites all modes stochastically. Hence, GWs can 
be detected if I = 2 mode amplitudes show a 
significantly different behavior than it would be 
expected from excitation by turbulent convection 
only; in particular, the ideal case is dominance of 
GW excitation, that is, I ^ 2 mode amplitudes are 
surmounted by amplitudes of / = 2 modes. If this 
is not the case, one can deduce an upper bound 



on an SBGW by employing Equations (64) and 



( 65 1 and assuming that excitation of I = 2 oscilla- 



tion modes is entirely due to excitation by GWs. 
Even if it turns out that the former case (the Sun 
as a detector for GWs) does not apply, the latter 
case would still be very interesting, given the fact 
that to date there is no method and experiment 
to derive a direct and significant upper bound on 
an SBGW at intermediate frequencies between the 
pulsar bound at low frequencies around 10 -8 Hz 
and the bounds by Earth-based interferometric de- 
tectors at high frequencies around 100 Hz. This in- 
termediate frequency region would be particularly 
interesting in the context of cosmic string models 
for an SBGW, since such a bound would signifi- 
cantly constrain possible values for e in small loop 
models with high values of J^gw for the constant 
part of the spectrum and therefore it would in par- 
ticular constrain the parameter space of small loop 
models, i.e., models where the size of string loops 
is determined by gravitational backreaction. 

In future work, we will apply the method de- 
veloped in this paper to stars other than the Sun 
focusing in particular on those stars, for which, 
due to a different stellar structure, the excita- 
tion mechanism proposed in the present paper 
might be even more effective compared to exci- 
tation by turbulent convection. For stars more 



massive than the Sun and having larger radii R, 
oscillation amplitudes will benefit from the effec- 
tive i? 2 -dependcnce present in the formulae for 
mean-square amplitudes and square surface veloc- 
ities (cf. Equations ((35), ([38]), and Q). How- 
ever, making definite further statements about 
other stars in this regard is intricate, since de- 
tailed knowledge about the eigenfunctions of a star 
is needed to calculate the resonance factors, which 
determine the strength of mean-square amplitudes 
and rms surface velocities. In combination with 
experimental data from the Kepler mission, the 
application of the method developed in this paper 
to stars other than the Sun might enable us to ei- 
ther use stars as giant hydrodynamic detectors for 
GWs, and this would set the scene for the detec- 
tion of GWs by asteroseismology, or it might en- 
able us to deduce a significant direct upper bound 
on an SBGW at intermediate frequencies. 

The authors are grateful to M.-A. Dupret and 
K. Belkacem for sharing their g-mode damping 



rates published in Belkacem et al. (2009) and 
to J. Christensen-Dalsgaard for the provision of 
the Aarhus adiabatic pulsation package ADIPLS. 
The authors also wish to thank O. von der Liihe 
for valuable discussions. M.R. acknowledges sup- 
port from the European Helio- and Asteroseismol- 
ogy Network (HELAS), which was funded as Co- 
ordination Action under the European Commis- 
sion's Framework Programme 6. D.M.S. acknowl- 
edges support from the German National Aca- 
demic Foundation. 



20 



A. Derivation of the Euler equations 



This appendix is concerned with explicitly deriving the equation of continuity ( 10 1 and the equation of 



motion ( |ll| for the internal velocity field from the general relativistic conservation equation of energy and 
momentum ([2]). 

As described in Sections 2.1 and 2.2 the following assumptions and approximations are in order: 

(i) Long-wavelength approximation. The GW tensor and hence by Equation ([8]) the Riemann tensor 
have no spatial dependence over the volume of the star. 



(ii) TT gauge. The GW tensor can be assumed to be in TT gauge (cf. Equation Q) when computing 
components of the Riemann tensor due to the invariance property of the Riemann tensor. 

(hi) Newtonian internal motions. Components of the internal velocity field v of the star are small compared 
to the speed of light, v % <C c. 

(iv) Linearization. Only terms linear in /i M „ need be retained (linearized theory assumption); in addition, 
only terms linear in the coordinates are retained. 



Furthermore, we define the following two tensors for convenience: 

( u^u v Y 
T» v = pu»u" + p — =- + g» v = M» v + P^. 



(Al) 



The components of the four-velocity field are given by it M = 7(0, v), where 7 = l/yT+ v 2 /c 2 . 
A.l. Equation of continuity 



The equation of continuity ( 10 ) follows from the zero-component of the general relativistic conservation 
equation of energy and momentum, 

d^M^ + = -T° 0Q T 00 - r^T 00 , (A2) 

where we only retained terms including T 00 on the right-hand side due to (hi). 
Using 7 ~ 1 in accordance with (iii) , we calculate 

= 9 (7pc 2 )+9 i ( 7 Vc) 



dp 
dt 



V-(pv) 



(A3) 



Furthermore, according to (i)-(iv) and Equations we have 

= (d p)( 1 2 + g 00 ) + (d lP )( 1 



+p 



d j 2 + di 7 



c 

+ V° 



0, 



since g 



on 



'1, 



0, and 



R°\) x k = 0. 



(A4) 



(A5) 



The first term in brackets vanishes, i? n fc ! ' i = 0, as can be immediately seen by employing Equation (|8j). The 
second term in brackets is also readily seen to be zero by using Equation m\ and applying (ii), R°^ k = 0. 
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The Christoffel symbols 



(A6) 



are evaluated in Fermi normal coordinates by using expressions Q-Q. The symbols needed in Equation 
(A2) are easily obtained as follows: 



1 oo 



^5°°<9 .9oo + 9 m {d g t - ^500) - 0. 



This is due to g 0z and dogoo already being of second order in the coordinates: g° l ~ and 



Furthermore, 



<9o5oo - -Rokoidoix 'x ) = 



= ^g 00 d goo + g 0t d goi + \g %0 d o g t3 
~ 0, 



(A7) 
(A8) 



where we additionally used 



dogij 



-R ikjl d {x k x l ) = 0. 



Inserting expressions (A3), (A4|, (A7), and (A9| into Equation (A2) finally yields the desired result, 



dp 

at 



V-(pv) = 0, 



(A9) 
(A10) 

(All) 



which is the usual equation of continuity. 

A. 2. The Euler equation 

The equation of motion for the internal velocity field is obtained from the spatial components of the 
general relativistic conservation equation of energy and momentum, 



d^M"* + d^p* = -r 00 T 00 , 

where, again, only T 00 terms were retained on the right-hand side due to (iii) 
As before, we start by evaluating the quantities on the left-hand side: 

a ( ( o7 2 CT 4 ) + a J (p7 2 ^^') 



(A12) 



dv % 



dp 

dt 



dv % 
~dt 



(v ■ vy 



(A13) 



where (iii) and Equation (All) were employed. By making use of (i), (iii), (iv), and Equations Q-^, one 
easily finds for the pressure tensor 



V" = d p 7 A -+g Ut )+d kP 7 



+ d k 7 



d (j 2 - 
ifd k p + pd^ 



-g 



(A14) 
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where in the second step we solely used (iii) and (iv). For the third step, the following expressions are 
needed: R\ l ■ — and R 3 j k — 0. The former can be read off Equation QSj) and the latter is obtained from 



Equation (18]) by applying (i). Hence, it is easily seen that 



0, 



(A15) 



and thus Equation (A14) follows 



For the right-hand side of Equation ( A12 1 one needs to evaluate 



1 



r oo = 9 la do9o* - ;A.9oo = ~^9 13d j9oo 



1 



1 



-—h* . x j 
2c 2 J ' 



where in the second step we made use of Equation ( A8 ) and 



d Q g z - -^Roktid Q (x k x l ) = 0. 



(A16) 



(A17) 



In the fourth step of Equation (A16), Rawi = — ^ih^i has to be inserted, which is readily obtained from 
Equation Q by assuming (ii). 



Thus, by inserting expressions (|A13|), (|A14|), and (|A16| into Equation (A12) we finally arrive at 

-Vp + / GW , 



P 



c)v 



at 



(v • V)v 



where the additional force per unit volume / GW is defined by the components 

/gw = nP^jX 3 ■ 



(A18) 



(A19) 



Employing Equation (All) and using tensor product notation, Equation (A18) is written as 

dpv 



dt 



V-(pv ® v) = -Vp + / Gw , 



(A20) 



which is the equation of motion to be derived in this subsection. 



B. Some properties of the gravitational wave tensor in transverse-traceless gauge 



In order to clearly state the relevant definitions and equations employed in Sections |3.1| and |3.2| we 
dedicate this appendix to a very brief review of some of the basic properties of the GW tensor in TT 
gauge. For a textbook presentation of the subject see, e.g., Maggiore (2008). 



In TT gauge (cf. Equation (19j)), the GW tensor satisfies the linearized field equations in vacuo of the form 



Dhij = 0. 



(Bl) 
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Since hij is traceless and symmetric (spin-2 operator), it can be expanded as 



2m 

ij > 



(B2) 



where the expansion coefficients h m are called the five independent spherical components of hij and where 
the expansion was made in terms of the following basis for the five-dimensional space of traceless symmetric 
tensors: 



y 



2, ±2 



y 



2,±1 




1 


±i 




±i 


-1 


i) 















1 


(• 





±i 




±i 





-1 





\ 





-1 











2 / 



= y?Pn i n j , 



y 1 



These tensors are related to spherical harmonics by the following relation: 

^2m( 

where n, denote the components of the radial unit vector 

n = e r = (sin 9 cos (f>, sin 8 sin cj>, cos 9) . 

Furthermore, we mention the orthogonality relation that applies to these tensors: 

15 



1:1 



With the help of this orthogonality relation, Equation (B2| can be inverted to give 

8tt. 



15 



ha (y?p) 



(B3) 
(B4) 
(B5) 

(B6) 
(B7) 

(B8) 
(B9) 



Finally, we discuss the plane wave expansion of the GW tensor in TT gauge. Due to Equation (Bl), hij 
can be expanded as 

hij(x) = J d 3 k [A tJ (ky kx + A* 3 {k)e- lkx ] , (BIO) 

where x = (ct, x) is used to denote spacetime points and where = (w/c, k) is the usual four- wave vector. 
The integral in Equation (BIO) can alternatively be expressed as an integral over frequency, v — lu/2tt, and 
solid angle, f2, 

1 



hij(x) 



dv v 2 I dQ 



c.c. 



(Bll) 



since d 3 k = \k\ 2 d\k\dil — (2-k j c) z v 2 dvd^L. Here n was used to denote the propagation direction of the 
individual plane waves, n = k/\k\. We now define the polarization tensors 



4(«) 

efjin) 



(B12) 
(B13) 
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where n given by Equation (B7), it, and v are pairwise orthonormal vectors, 



u = ee = (cos 9 cos <j), cos sin </>, — sin O) 
v = = (— sin <fi, cos <f>, 0). 



(B14) 
(B15) 



From this definition, it is clear that the polarization tensors represent transversal polarization states and 
that these polarization states are independent. The latter point is also reflected by the fact that 



zf i (n)e A '' i i(n) = 25 A ' A ', 



(B16) 



which can be proved explic itly. Here A = {+, x} labels the polarization states. Since applying the gauge 
conditions (|9j) to Equation (BIO I yields k l Aij(k) = 0, the tensor components Aij can be expanded in terms 
of the polarization tensors, 

;/ 2 



(B17) 



where the expansion coefficients fiA^, n) are called the Fourier amplitudes of hij. Substituting Equation 
(B17) into Equation (Bll) and extending the frequency integration domain to minus infinity by setting 
hA(—v.n) = h* A (v, n) finally yields 



hjj (t, x) = / dv I dQliA(v, n)e^(ri)e 

A J — OO J 



i27rv(t—n-x / c) 



(B18) 



This is the desired form of the plane wave expansion for the GW tensor employed in Section |3.2| 
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